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Abstract 

Introduction: Epithelial-mesenchymal transition (EMT) and nnesenchynnal-epithelial transition (MET) facilitate breast 
cancer (BC) metastasis; however, stable molecular changes that result as a consequence of these processes remain 
poorly defined. Therefore, with the hope of targeting unique aspects of metastatic tumor outgrowth, we sought to 
identify molecular markers that could identify tumor cells that had completed the EMTMET cycle. 

Methods: An in vivo reporter system for epithelial cadherin (E-cad) expression was used to quantify its regulation in 
metastatic BC cells during primary and metastatic tumor growth. Exogenous addition of transforming growth factor 
(31 (TGF-|3l) was used to induce EMT in an in situ model of BC. Microarray analysis was employed to examine gene 
expression changes in cells chronically treated with and withdrawn from TGF-|3l, thus completing one full EMTMET 
cycle. Changes in fibroblast growth factor receptor type 1 (FGFRl) isoform expression were validated using PCR 
analyses of patient-derived tumor tissues versus matched normal tissues. FGFRl gene expression was manipulated 
using short hairpin RNA depletion and cDNA rescue. Preclinical pharmacological inhibition of FGFR kinase was 
employed using the orally available compound BGJ-398. 

Results: Metastatic BC cells undergo spontaneous downregulation of E-cad during primary tumor growth, and its 
expression subsequently returns following initiation of metastatic outgrowth. Exogenous exposure to TGF-|3l 
was sufficient to drive the metastasis of an otherwise in situ model of BC and was similarly associated with a depletion 
and return of E-cad expression during metastatic progression. BC cells treated and withdrawn from TGF-|3 stably 
upregulate a truncated FGFR1-|3 splice variant that lacks the outermost extracellular immunoglobulin domain. 
Identification of this FGFRl splice variant was verified in metastatic human BC cell lines and patient-derived 
tumor samples. Expression of FGFRl -(3 was also dominant in a model of metastatic outgrowth where depletion 
of FGFRl and pharmacologic inhibition of FGFR kinase activity both inhibited pulmonary tumor outgrowth. 
Highlighting the dichotomous nature of FGFR splice variants and recombinant expression of full-length FGFRl-a 
also blocked pulmonary tumor outgrowth. 

Conclusion: The results of our study strongly suggest that FGFRl -(3 is required for the pulmonary outgrowth of 
metastatic BC. Moreover, FGFRl isoform expression can be used as a predictive biomarker for therapeutic application 
of its kinase inhibitors. 
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Introduction 

The reported results from several recent studies suggest 
that metastatic breast cancer (BC) cells undergo epithe- 
lial-mesenchymal transition (EMT) during invasion and 
dissemination and that the reverse process of mesenchy- 
mal-epithelial transition (MET) occurs at some point dur- 
ing metastatic tumor outgrowth [1-3]. In fact, the ability 
of BCs to transition between an epithelial and mesenchy- 
mal state seems to be a key feature of the metastatic 
process and has recently been more accurately termed epi- 
thelial-mesenchymal plasticity [4]. Indeed, many of the 
well-established changes in gene expression that take 
place during EMT return to baseline during MET. How- 
ever, stable phenotypic markers capable of distinguishing 
metastatic BC cells that have undergone an EMTMET 
cycle from their indolent counterparts that have not 
undergone an EMT program remain to be identified. 

Recently, we demonstrated that transforming growth 
factor p (TGF-p)-induced EMT empowers BCs with the 
ability to invade in response to paracrine epidermal 
growth factor (EGF) stimulation, thereby facilitating the 
egress of BC cells from the primary tumor [5]. Intri- 
guingly, we also observed that cells proficient in under- 
going metastatic outgrowth downregulate EGF receptor 
(EGFR) during the MET process [2]. These findings may 
serve to explain the contrasting clinical data that estab- 
lish EGFR as a predictor of poor prognosis for BC pa- 
tients [6], yet administration of monotherapies directed 
against EGFR or in conjunction with other chemother- 
apies has failed to provide a clinical benefit for BC pa- 
tients [7-9]. 

Increased expression of fibroblast growth factor recep- 
tor (FGFR) types 1 and 3 have recently been identified 
as two of six receptor tyrosine kinases associated with 
poor disease-free survival and/or decreased overall sur- 
vival in BC patients [10]. FGFRl, FGFR2 and FGFR3 all 
exist as several different isoforms generated via alterna- 
tive splicing [11]. Two of the best-described variants are 
generated via inclusion of unique versions of the third 
(III) extracellular immunoglobulin (Ill-Ig) domain, and 
hence they are termed FGFR-IIIb and FGFR-IIIc [11]. 
The Ill-Ig domain governs the specificity of FGFR bind- 
ing to the 18 different FGF ligands. FGF2 (basic FGF), 
for instance, has an extremely high affinity for the IIIc 
isoform [12]. Another FGFR splicing event results in ei- 
ther the inclusion (FGFR-a) or exclusion (FGFR-p) of 
the first Ig domain and/or the linker region between Igl 
and Igll, an area termed the acid box [13]. Importantly, 
Igl and the linker region regulate the affinity of FGFR 
for its particular ligand [14]. Furthermore, increased ex- 
pression of the p versus a isoform of FGFRl has been 
correlated with reduced relapse-free survival in a cohort 
of BC patients [15]. The use of antisense morpholino oli- 
gonucleotides to convert FGFR splicing from the p to 



the a isoform can induce apoptosis in glioblastomas 
[16]. Currently, the upstream mechanisms that regulate 
FGFRa/FGFRp splicing remain poorly defined, but TGF- 
p and its induction of EMT can cause upregulation of 
FGFRl-IIIc and downregulation of FGFR2-IIIb [17]. How- 
ever, little is known about the mechanism by which either 
of these events drives the metastatic progression of BCs. 

In The present study, we sought to identify factors that 
are stably altered during EMT:MET cycles and thus 
might act as drivers of metastatic tumor outgrowth. To 
this end, we utilized microarray expression analyses of 
BC cells that had been treated and withdrawn from ex- 
ogenous TGF-p to uncover a stable upregulation of the 
FGFRl -p isoform, an event that was also readily detected 
in samples obtained from BC patients and in metastatic 
human BC cell lines. Along these lines, genetic depletion 
of total FGFRl and/or ectopic overexpression of FGFRl- 
a in the D2.A1 model of pulmonary metastatic out- 
growth potently inhibited pulmonary tumor formation. 
Collectively, these studies establish FGFRl -p as a critical 
player whose expression is stably altered in metastatic 
BCs that have experienced oncogenic TGF-p signaling 
and undergone EMTiMET cycles. Our findings highlight 
the need to further elucidate the pro- and antitumori- 
genic nature of FGFR to appropriately administer small- 
molecule inhibitors for the treatment of metastatic BC. 

Methods 

Cell lines and reagents 

Murine D2.A1 and human MCF-lOA derivatives (MCF- 
lOA, TIK and Calh) were ethically obtained from Dr 
Fred Miller (Wayne State University, Detroit, MI, USA) 
[18,19], and murine 4T1, human MCF-7 and human 
MDA-MB-231 cells were purchased from the American 
Type Culture Collection (ATCC; Manassas, VA, USA). 
In addition, the following human BC cell lines were ori- 
ginally purchased from the ATCC and ethically obtained 
from Dr John J Pink (Case Western Reserve University, 
Cleveland, OH, USA): MDA-MB-361, ZR-75-1, T47D, 
T47D-C42W and BT549 [20]. The aforementioned hu- 
man and murine BC cells were cultured in Dulbeccos 
modified Eagles medium supplemented with 10% fetal 
bovine serum and 1% penicillin-streptomycin as de- 
scribed previously [5]. Bioluminescent 4T1 and D2.A1 
cells were engineered to stably express firefly luciferase 
under the selection of Zeocin (InvivoGen, San Diego, 
CA, USA) as described previously [2,5,21]. Biolumines- 
cent normal mammary epithelial (NME) cells were 
constructed and cultured as previously described [5]. 
Dual-bioluminescent E-cad reporter cells were gener- 
ated by stably transfecting 4T1 cells with pcDNA3.1/ 
Hygro mammalian expression vector (Invitrogen, Carlsbad, 
CA, USA) encoding Renilla luciferase under control of the 
cytomegalovirus (CMV) promoter and pGL4.20[/Mc2/Puro] 
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vector (Promega, Madison, WI, USA) that encodes firefly 
luciferase under control of the human Cdhl promoter [22]. 
Cellular depletion of FGFRl expression was achieved by 
glycoprotein of vesicular stomatitis virus lentiviral trans- 
duction of TRC pLKO.l short hairpin RNA (shRNA) vec- 
tors (Thermo Scientific, Pittsburgh, PA, USA) (Additional 
file 1: Table SI) as described previously [2,21]. Ectopic ex- 
pression of FGFRl-a-IIIc was accomplished as described 
previously and selected for using neomycin [2]. 

In vivo bioluminescence imaging of tumor growth and 
metastasis 

Parental (that is, scrambled shRNA) and FGFRl- 
manipulated D2.A1 cells were injected into the lateral 
tail veins of 5-week-old female BALB/C mice (The 
Jackson Laboratory, Bar Harbor, ME, USA). Where 
indicated, mice bearing D2.A1 pulmonary tumors were 
treated daily with BGJ-398 (ChemieTek, Indianapolis, IN, 
USA) or PF-573271 (PF271; Pfizer Pharmaceuticals, New 
York, NY, USA) at 50 mg/kg by oral gavage. Alternatively, 
Cdhl reporter 4T1 cells (1 x 10^ cells) were engrafted onto 
the mammary fat pads of 4-week-old BALB/c mice. Circu- 
lating 4T1 tumor cells were isolated from the inferior vena 
cava at the time of necropsy using 3% sodium citrate. Fol- 
lowing lysis of red blood cells, circulating tumor cells were 
selected for with 5 (ig/ml Zeocin (the selectable marker 
for firefly luciferase). Luciferase-expressing NME cells 
(1 to 2 X 10^ cells) were engrafted onto the mammary 
fat pads of 5-week-old female nulnu mice. All bio- 
luminescent images were captured using a Xenogen 
IVIS 200 preclinical imaging system (Caliper Life Sciences/ 
PerkinElmer, Hopkinton, MA, USA) within the Small Ani- 
mal Imaging Resource Center at the Case Comprehensive 
Cancer Center as previously described [5,21,23]. 

Gene expression profiling 

NME cells were cultured in the presence of TGF-|3l 
(5 ng/ml) for 4 weeks, at the end of which TGF-pl sup- 
plementation was discontinued and the cells were allowed 
to recover for an additional 4 weeks. Total RNA was pre- 
pared from unstimulated cells of similar passage (pre- 
TGF) and the post-TGF NME ceUs. Microarray analyses 
were performed in triplicate using the GeneChip 
Mouse Gene ST 1.0 Array (Affymetrix, Santa Clara, 
CA, USA). Genes regulated more than twofold are 
given in Additional file 2: Table S2. The complete data 
set has been deposited in the National Center for Biotech- 
nology Information Gene Expression Omnibus (GEO) 
database [GEO:GSE54491] [24]. 

mRNA transcript analyses 

For real-time PCR analysis, normal murine mammary 
gland (NMuMG) and NME ceUs were stimulated with 
TGF-|3l (5 ng/ml) for varying lengths of time, and then 



total RNA was isolated using RNeasy Plus Mini Kit 
(QIAGEN, Valencia, CA, USA). Afterward, total RNA 
was reverse-transcribed using the iScript cDNA Synthe- 
sis Kit (Bio-Rad Laboratories, Hercules, CA, USA), and 
semiquantitative real-time PCR was conducted using 
iQ SYBR Green Supermix (Bio-Rad Laboratories) as 
described previously [21]. Identification of FGFR spUce 
variants was accomplished by visualizing PCR products 
separated by gel electrophoresis. The oligonucleotide 
primer pairs used are provided in Additional file 1: 
Table SI. 

Immunoblotting and immunohistochemical analyses 

For immunoblot assays, lysates generated from two- and 
three-dimensional cultures were prepared as described 
previously [2]. The antibodies we used are described in 
Additional file 3: Table S3. For immunohistochemistry, 
tissues were fixed in 10% formalin, and histological sec- 
tions were prepared by the Tissue Procurement and 
Histology Core at the Case Comprehensive Cancer Cen- 
ter. Sections were deparaffinized and stained with the in- 
dicated antibodies given in Additional file 3: Table S3. 

Three-dimensional organotypic growth assays 

Cells were diluted in complete media supplemented with 
5% Cultrex reagent (Trevigen, Gaithersburg, MD, USA) 
and seeded onto solidified Cultrex cushions (50 (il/well) 
contained in 96-well plates (1 x 10^ cells/cm^). Longitu- 
dinal bioluminescence growth assays were performed as 
described previously [2,5]. Pharmacological inhibitors 
against FGFR (BGJ-398) or focal adhesion kinase (FAK) 
(PF271) were added to cultures at the indicated concen- 
trations and times. 

Study approval 

All animal procedures were performed in accordance 
with protocols approved by the Institutional Animal 
Care and Use Committee of the Case Western Reserve 
University School of Medicine. For human BC speci- 
mens, primary tumors and matched normal tissues were 
collected and processed under protocols approved by 
the Institutional Review Board of the Cleveland Clinic. 
All patients provided their written informed consent 
allowing the study investigators to have access to their 
tumor specimens and clinical data. 

Statistical analyses 

Statistical values were defined using an unpaired Stu- 
dents ^-test, where a P-value less than 0.05 was consid- 
ered significant. Statistically significant differences in the 
overall survival of mice bearing control and FGFR- 
manipulated D2.A1 pulmonary tumors were analyzed 
using a logrank test. P values for all experiments are 
indicated. 
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Results 

E-cadherin is dynamically regulated during spontaneous 
breast cancer metastasis 

We recently established that BCs that can downregulate 
E-cadherin (E-cad) are at a selective advantage to initiate 
metastatic outgrowth within the pulmonary microenvir- 
onment [2,25]. However, these and other studies suggest 
that, as pulmonary metastases progress from micro- to 
macroscopic lesions, E-cad expression returns in a man- 
ner reminiscent of the differentiation status observed in 
noninvasive primary tumors [1,2]. Along these lines, the 
highly metastatic 4T1 model of late-stage BC displays 
several epithelial characteristics, including the expression 
of E-cad [2,21]. Therefore, we sought to use this model 
system in combination with dual-substrate biolumines- 
cence [26] to track the in vivo dynamics of E-cad regu- 
lation as 4T1 cells progress through the metastatic 



cascade. To this end, 4T1 cells that stably expressed 
Renilla luciferase under the control of the CMV pro- 
moter were engineered to stably express firefly lucifer- 
ase under the control of the E-cad promoter [2]. Using 
this approach, we observed areas within primary tumors 
that failed to express E-cad luminescence (Figure lA). 
Large areas in the center of the primary tumor were nec- 
rotic and therefore lacked luminescence, so we utilized 
immunohistochemistry (IHC) to verify a lack of E-cad 
protein in several outer and viable areas of the tumor 
(Figure lA). Furthermore, early nodular metastases 
produced little to no firefly luminescence, and they 
were similarly verified by IHC analysis to be viable 
tumor tissue and E-cad-negative (Figure lA). Once 
these metastatic lesions had spread deeper into the 
pulmonary tissue, however, E-cad luminescence and 
protein returned in a manner consistent with that of a 
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Figure 1 Spontaneously metastatic breast cancer cells dynamically regulate epithelial cadherin during in vivo primary and metastatic 
tumor growth. (A) Bioluminescence imaging (BLI) of metastatic 4T1 cells expressing cytomegalovirus (CM\/)-driven Renilla luciferase and firefly 
luciferase driven by the epithelial cadherin (E-cad) promoter. Cells were engrafted onto the mammary fat pads of BALB/c mice, and qualitative 
E-cad promoter activity (for example, luciferin-derived bioluminescence) was monitored in vivo and ex vivo. E-cad bioluminescence was spatially 
correlated with E-cad protein expression as determined by IHC. Bars indicate 40x, lOOx and 400x magnifications. (B) E-cad promoter activity (for example, 
luciferin-derived bioluminescence) was quantified by normalization to CMV promoter activity (for example, coelenterazine-derived bioluminescence) in in 
vitro cultured cells and ex vivo tissues derived from primary tumors and their late-stage metastases (Mets). Data represent ten samples (n = 10, ±SE) 
derived from five individual mice bearing primary and metastatic tumors, resulting in the indicated lvalues. (C) Circulating tumor cells (CTCs) were 
isolated from the blood of 4T1 tumor-bearing mice and selected for resistance to Zeocin. Photomicrographs show the antibiotic-resistant CTCs (CTC 
exwVo). The morphologically epithelial cells (black arrow) were physically isolated from the mesenchymal cells (white arrow) and further 
subcultured as separate populations of mesenchymal-like and epithelial-like cells as indicated. Bars indicate a lOOx magnification. 
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MET program (Figure lA). To verify that the differ- 
ences in firefly bioluminescence were due to regulation 
of the E-cad promoter, primary and metastatic tumors 
were homogenized and the firefly luminescence was 
normalized to the CMV-driven expression of Renilla 
luciferase expressed in the same cells. This quantitative 
approach resulted in an 8.59-fold upregulation of E-cad 
promoter activity in pulmonary metastases as compared 
to their corresponding primary tumors (Figure IB). Con- 
sistent with the role of EMT in facilitating primary tumor 
dissemination, we readily captured highly mesenchymal- 
appearing, antibiotic- resistant tumor cells from the per- 
ipheral blood of mice bearing primary 4T1 tumors 
(Figure IC). When left as a general cell population, 
however, the circulating tumor cells quickly resumed 
an epithelial phenotype in culture (Figure IC). Taken 
together, the results of these studies suggest that the 
ability to efficiently transit between epithelial and mes- 
enchymal states is associated with efficient completion 
of the metastatic cascade. 

TGF-P treatment Is sufficient to drive orthotopic mammary 
tumor metastasis 

TGF-p is one of several "master regulators" of the EMT 
process [27]. Given the spontaneous nature of E-cad 
regulation during 4T1 metastasis, we next sought to spe- 
cifically address whether TGF-p treatment would be suffi- 
cient to facilitate the metastasis of primary orthotopic 
mammary tumors. To do so, we utilized our established 
model in which NMuMG cells are transformed via over- 
expression of EGFR (NME cells) [5]. Although orthotopic 
tumor formation in this model is robust, postsurgical re- 
currence and distant metastasis were not observed 
(Figure 2A). In stark contrast to control NME tumors, 
however, treatment of NME cells with TGF-p prior to 
mammary fat pad engraftment (Figure 2B) yielded larger 
primary tumors (Figures 2A to 2D) that displayed in vivo 
modulation of the EMT markers E-cad and fibronectin, as 
well as downregulation of estrogen receptor a (ER-a) 
(Figure 2E). More importantly, these TGF-p-treated tu- 
mors displayed aggressive recurrence following primary 
tumor resection (Figures 2A, 2F and 2G), and progressed 
to form pulmonary metastases in a fashion that was sig- 
nificantly increased as compared to the control tumors 
(Figure 2H). Taken together, these data clearly establish 
that TGF-p treatment is sufficient to induce the postsurgi- 
cal recurrence and metastasis of an otherwise indolent 
model of in situ mammary carcinoma. 

FGFR is stably upregulated following TGF-P treatment 
and withdrawal 

Given the transient nature of exogenous TGF-p treat- 
ment prior to NME cell engraftment, we next sought to 
assess the regulation of MET in this model. Consistent 



with our IHC analysis (Figure 2E), ex vivo subculture of 
post-EMT NME primary tumors displayed a highly mes- 
enchymal phenotype in two- and three-dimensional cul- 
tures and downregulation of E-cad (Additional file 4: 
Figures SIA and SIB). TGF-p-treated NME tumor cells 
that had undergone metastasis were subcultured from 
the lungs of mice. In comparison to their parental NME 
counterparts, these cells underwent dramatically enhanced 
primary tumor formation, postsurgical recurrence and 
spontaneous pulmonary metastasis upon secondary mam- 
mary fat pad engraftment (Additional file 4: Figure SIC). 
These resulting metastases were subcultured and termed 
the NME lung metastatic (NME-LM2) cell line (Additional 
file 5: Figure S2A). Further in vitro analyses of the NME- 
LM cell lines revealed a return of E-cad expression to levels 
that approximated those detected in their untreated NME 
parental cells (Additional file 5: Figures S2B and S2C). 
Taken together, these results indicate that depletion of 
E-cad by exogenous TGF-p treatment was stably main- 
tained during formation of the primary tumor but 
readily returned to baseline expression levels during 
formation of macroscopic pulmonary metastases. 

Given that TGF-p treatment and metastasis are cap- 
able of selecting for a spontaneously metastatic cell 
model that has renewed E-cad expression, we next sought 
to identify^ stable changes in gene expression that could 
characterize cells that reside in either a pre- or post-TGF- 
p exposure state. To do so, we conducted microarray 
analyses of NME cells that had undergone 4 weeks of 
exogenous TGF-pl treatment, which was followed by 
an additional 4- week withdrawal of exogenous TGF-pl 
(Additional file 2: Table S2). Indeed, following this ex- 
perimental protocol, the expression patterns of trad- 
itional EMT markers, such as E-cad, N-cadherin, vimentin, 
fibronectin, a-smooth muscle actin. Twist and Snail all 
returned to baseline levels (Additional file 2: Table S2). 
Surprisingly, 98 genes were stably modulated more 
than threefold following this TGF-p treatment and 
withdrawal protocol. Gene set enrichment analysis re- 
vealed that the highest degree of overlap between 
genes upregulated in our gene list was that of genes 
downregulated in mouse embryonic fibroblasts after 
10 hours of TGF-p treatment [GEO:GSE15871] [28]. 
Furthermore, 15 upregulated and 28 downregulated 
genes in our data set were shared with a data set gen- 
erated upon knockout of the TGF-P family member 
bone morphogenic protein 2 [29]. The transcription 
factor Snai2, which functions to inhibit the expression 
of E-cad [30], was downregulated tenfold in cells that 
had undergone this cycle of TGF-p treatment and 
withdrawal. Taken together, these findings suggest that 
certain aspects of TGF-P signaling can be stably al- 
tered, leading to a unique transcriptional profile in cells 
following prolonged ligand exposure and withdrawal. 
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Figure 2 Transforming growthi factor P treatment is sufficient to induce postsurgical recurrence and pulmonary metastasis. (A) Normal 
mammary epithelial (NME) cells were left untreated (before transforming growth factor treatment (pre-TGF)) or treated with TGF-(3l (post-TGF) 
and engrafted (2x10^ cells/mouse) onto the mammary fat pad, at which point primary tumor growth was monitored by bioluminescence. Five 
weeks after engraftment of pre- and post-TGF NME cells, the primary tumors were surgically removed and tumor recurrence was similarly monitored 
by bioluminescence imaging. Shown are two representative mice from each experimental cohort, both before and 3 weeks after surgical resection of 
the primary tumor. (B) NME cells were treated with TGF-(3l (Post) as described in (A) and monitored for expression of the mesenchymal marker 
fibronectin (FN) and the marker epithelial cadherin (E-cad). Estrogen receptor a (ERa) was potently downregulated in NME cells following 
TGF-(3 treatment. (C) Photograph showing ex vivo pre- and post-TGF primary tumors following surgical resection demonstrating 
complete and intact primary tumor excision. (D) The average weight (±SE) of pre- and post-TGF primar/ tumors shown in (C) and the indicated 
P value. (E) Pre- and post-TGF primary tumors were analyzed by immunohistochemistry for the in vitro markers described in Panel B. All images were 
taken at lOOx magnification and insets were taken at 400x magnification. (F) Growth of pre- and post-TGF NME primary and recurrent tumors were 
monitored using digital caliper measurements (n=5 mice/group). Vertical line indicates time of primary tumor resection. Data are the mean (±SE) 
tumor size for each group . *P < 0.05. (G) Photograph of a representative mouse 3 weeks after resection of a post-TGF NME tumor. Yellow 
outlines show widespread recurrent tumor formation. As shown in (A), mice bearing pre-TGF NME tumors failed to form recurrent tumors. 
(H) Pulmonary metastasis of pre- and post-TGF NME tumors was quantified by bioluminescence. Line indicates time of primary tumor 
resection. Inset: Bioluminescent images of ex vivo lungs from representative mice bearing pre- and post-TGF NME tumors. Data are the mean 
(±SE) area flux values for each group. ""P < 0.05 (n = 5 mice/group). 



Further analysis of the microarray data established that 
three EGFR ligands, betacellulin, amphiregulin and epir- 
egulin were all potently downregulated (Additional file 6: 
Figure S3). Conversely, FGFRl, FGFR2 and FGFR3 



were upregulated 15. 8-, 4.2- and 14.8-fold, respectively 
(Figure 3A). Given these findings and a growing interest 
in targeting FGFR as a treatment for BC [31], we sought 
to further assess the role of FGFR in facilitating 
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Figure 3 Upregulation of fibroblast growtli factor receptor type 1 is stably maintained during oncogenic transforming growth factor (3 
signaling. (A) Normal mammary epithelial (NME) cells were left untreated (before transforming growth factor treatment (pre-TGF)) or treated and 
allowed to recover from exogenous TGF-(3l treatment (post-TGF-Rec) as described in the Methods section. Global gene expression was assessed 
by microarray analyses. Data are results of RT-PCR analysis confirming stable upregulation of fibroblast growth factors 1 through 3 (FGFRl to 
FGFR3) following TGF-(3l treatment and withdrawal. (B) Real-time PGR analyses showing the differential expression of FGFRl to FGFR4 in control 
yellow fluorescent protein (YFP), nonmetastatic (NME) and lung metastatic (LM2) cell lines derived from the normal murine mammary gland (NMuMG) 
cells. (C) Control NMuMG (YFP) and NME cells were not stimulated (NS), stimulated with TGF-(3l (5 ng/ml) for 48 hours (JGF) and subsequently recovered 
for an additional 48 hours without exogenous ligand (Rec). Expression of FGFRl was analyzed by real-time RT-PCR. Data in panels B and C are the mean 
expression values (±SD) of three independent experiments, resulting in the indicated P-values. (D) Human MCF-1 OA-derived breast cancer (BC) cells 
of increasing metastatic potential were analyzed by RT-PCR for their expression of FGFRl . Glyceraldehyde 3-phosphate dehydrogenase 
(GAPDH) served as a loading control. Data are representative of three independent experiments. (E) In silico analyses of the Gene Expression 
Omnibus data set [GEO:GSE20437] demonstrating significant increases in FGFRl to FGFR4 in prophylactically removed breast tissue (Pro), as 
well as in ER-negative (ER-) and ER-positive (ER+) tumor samples, as compared to a cohort of normal (Nor) breast tissue samples. Data are 
the individual values for each sample resulting in the indicated mean values (+SE) and P values. 



pulmonary metastatic outgrowth. In doing so, we first ex- 
amined FGFR expression in our NME progression series. 
In this system generated in vivo, FGFRl was the predom- 
inant family member stably upregulated in the meta- 
static NME-LM2 cells (Figure 3B). Using a condensed 
in vitro TGF-p treatment and withdrawal experiment, 
we demonstrate the stable upregulation of FGFRl in 
the transformed NME cells as compared to their non- 
transformed NMuMG counterparts (Figure 3C). Simi- 
larly to the NME progression series, the human MCFIOA 
BC progression series also displayed a robust increase in 
FGFRl expression in cells with a higher metastatic poten- 
tial (Figure 3D). Importantly, analysis of 20 human BC 
samples revealed a five- to tenfold increase in the expres- 
sion of FGFRl to FGFR4 in tumor tissue as compared to 
matched normal tissue (Additional file 7: Table S4). 



Furthermore, in silico analysis of the GEO data set [GEO: 
GSE20437] in which we compared normal and prophylac- 
tically removed mammary tissue to ER + and ER - tumor 
tissues demonstrated significant increases in several FGFR 
family members (Figure 3E). In light of other recent find- 
ings that suggest a role for FGFRl in driving the progres- 
sion of basal-type BCs [10,32-34], our findings clearly 
indicate that FGFRl is upregulated in BCs with increasing 
metastatic potential. They also suggest that FGFRl may 
be capable of driving the pulmonary outgrowth of meta- 
static BCs. 

FGFRl -p-lllc is selected for during metastatic progression 

TGF-p has recently been shown to be capable to modu- 
lating FGFR expression and alternative splicing of the 
FGFRl-IIIb/FGFRl-IIIc exon [17]. Indeed, examination 
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Figure 4 The p isoform of the fibroblast growth factor receptor type 1 is selected for in increasingly metastatic cells and can be 
readily identified in patient tumor samples. (A) Schematic representation tine FGFRl transcript depicting tine coding regions for tine 
immunoglobulin (Ig), transmembrane (TM), and kinase domains. The location of the unique lllb/lllc primer sets (arrows) and the flanking primer 
set detecting the inclusion or exclusion of the a exon (arrows) are also indicated. (B) Expression of FGFRl isoforms depicted in (A) were analyzed 
by RT-PCR in the nonmetastatic normal mammary epithelial (NME) cells and their lung metastatic (LM) and isogenic NME-LM2 counterparts 
before (ns) and after a 48-hour treatment with 5 ng/ml transforming growth factor (3 (JGF). (C) Luminal MCF-7, T47D, MDA-MB-361 and ZR-75-1 and 
basal MDA-MB-231 and BT549 human breast cancer cells were analyzed for inclusion or exclusion of the a exon of FGFRl. The estrogen-independent 
T47D-C42W cells were also analyzed. Data in (B) and (C) are representative of three independent experiments. (D) Four of the patient tumor samples 
(T; patients 1, 2, 7 and 13 (PI, P2, P7 and PI 3, respectively) that demonstrated upregulation of FGFRl as compared to their matched normal mammary 
tissues (N) (Additional file 7: Table S4) were further analyzed for expression of the a versus (3 isoforms of FGFRl . The fold increase in total FGFRl 
expression as determined by real-time PGR analysis that was normalized to glyceraldehyde 3-phosphate dehydrogenase (GAP) is indicated. 



of NME cells using isoform-specific RT-PCR primer sets 
(Figure 4A) indicated that TGF-|3 treatment upregulated 
the expression of FGFRl-IIIc (Figure 4B). Importantly, 
this modulated expression of FGFRl -IIIc was maintained 
throughout the metastatic process as NME-LM2 cells dis- 
played constitutive upregulation of FGFRl-IIIc (Figure 4B). 
However, FGFRl is subject to further alternative splicing 
that leads to the inclusion or exclusion of the third (a) 
exon that encodes the outermost Ig domain (Igl) of the re- 
ceptor [12] (Figure 4A). Importantly, increased expression 
of the |3 isoform of FGFRl has previously been associated 
with breast cancer [15]. Although inclusion of the a exon 
is dominant following in vitro TGF-|3 treatment of NME 
cells, it is the truncated |3 isoform that was selected for 
and enriched during the metastatic generation of the 
NME-LM2 cells (Figure 4B). Examination of several hu- 
man EC cell lines revealed that the luminal MCF-7, T47D, 
MDA-MB-361 and ZR-75-1 cells all express little to no 
FGFR1-|3 (Figure 4C). In contrast, the basal-like and meta- 
static BT549 and MDA-MB-231 human BC cells demon- 
strate robust expression of the truncated FGFRl -|3 isoform 
(Figure 4C). Interestingly, examination of a T47D subline, 
C42W, that was selected for by prolonged estrogen 
deprivation [20] displayed a loss of the FGFRl -a isoform 
and a slight gain in FGFR1-|3 (Figure 4C). Furthermore, 



two of the four patient tumor samples that displayed in- 
creased FGFRl expression as compared to their matched 
surrounding normal tissues (Additional file 7: Table S4) 
also demonstrated expression of the truncated FGFRl-p 
isoform (Figure 4D). Collectively, these results strongly 
suggest that expression of FGFRl -|3 is increased with meta- 
static progression. 

Genetic depletion of FGFRl prevents pulmonary tumor 
outgrowth 

Our data thus far demonstrate that FGFRl is stably up- 
regulated following oncogenic TGF-p signaling (Figure 3A) 
and its stimulation of pulmonary metastasis (Figure 4B). 
To validate these findings in the context of metastatic out- 
growth, we next took a genetic approach to specifically 
target FGFRl expression in the well-established D2.A1 
model of pulmonary metastatic outgrowth [2,35,36] . In ac- 
cordance with the NME progression series (Figure 4B), 
the pulmonary outgrowth proficient D2.A1 cells readily 
express FGFRl -p-IIIc (Additional file 8: Figure S4). Using 
pan-FGFRl primers (Figure 5A), as well as those that de- 
pict FGFRl-a and -p isoform expression (Figure 5B), we 
identified three unique shRNAs that elicit differential de- 
pletion of FGFRl. All of the FGFRl -targeting constructs 
inhibited the outgrowth of D2.A1 cells when cultured in a 
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Figure 5 Depletion of fibroblast growth factor receptor type 1 inhibits pulmonary tumor outgrowth. (A) Real-time PCR analysis of fibroblast 
growth factor receptor type 1 (FGFRl) in D2.A1 cells expressing a nontargeting control short hairpin RNA (shRNA; scram) or three unique shRNA 
sequences targeting FGFRl (sh#1 to sh#3). Data are the mean (+SD) of three independent experiments. (B) RT-PCR analysis of D2.A1 cells 
expressing the FGFRl shRNA targeting sequences as described in (A). Primers used flank the a exon of FGFRl and thus depict inclusion (a) or 
exclusion (|3) of this exon. GAPDH, Glyceraldehyde 3-phosphate dehydrogenase; sc, scrambled short hairpin RNA. (C) Three-dimensional 
outgrowth of the FGFRl -depleted D2.A1 cells described in (A) and (B) was quantified by bioluminescence over the course of 11 days (DO 
to D1 1). Data are the means (±SE) of two independent experiments completed in triplicate, resulting in the indicated lvalue. RLU, Relative light units. 
(D) D2.A1 cells expressing FGFRl shRNAs #2 or #3 were injected into the lateral tail vein of female BALB/c mice, and pulmonary tumor formation was 
quantified by bioluminescence. Data are the mean (+SE) area flux values at the indicated time points expressed as a percentage of the injected value 
(% of TO; P < 0.01; n = 5 mice/group). Insets: images of representative lungs 21 days after tail vein injection with the indicated cells. (E) Pulmonary 
tissues from mice injected with control (scram) and FGFRl -depleted (shFr1#2) D2.A1 cells were stained with hematoxylin and eosin (H&E) or analyzed 
by immunohistochemistry for phosphorylation of extracellular signal-regulated kinases 1 and 2 (pErk1/2), the proliferation marker Ki67 or FGFRl. H&E 
images were taken at a 1x magnification, while all immunohistochemistry images were taken at a 400x magnification. (F) RT-PCR analysis of D2.A1 cells 
expressing a scrambled shRNA (sc) or the #2 shRNA targeting murine FGFRl (sh). Depletion of endogenous FGFRl in these cells was rescued by 
expression of a nontargeted human full-length form of FGFRl -a-lllc (Fr1-a) or green fluorescent protein (GFP) as a control. The relative ratio of a to 
(3 FGFRl transcripts is shown. (G) Three-dimensional outgrowth of the FGFRl -depleted and rescued D2.A1 cells described in (C) was quantified by 
bioluminescence. Data are the mean (±SE) of three independent experiments completed in triplicate, resulting in the indicated P value. 
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three dimensional (3D) matrix that recapitulates the pul- 
monary microenvironment (Figure 5C) [2,35-37]. Import- 
antly, the extent to which D2.A1 organoid outgrowth was 
inhibited was proportional to the level of FGFRl depletion 
(Figures 5 A to 5C). Further in vivo examination of the two 
most effective FGFRl targeting sequences demonstrated 
that depletion of FGFRl potently blocked pulmonary 
tumor outgrowth (Additional file 9: Figure S5; Figure 5D). 
Immunohistochemical analysis of these pulmonary tumors 
indicated that depletion of FGFRl was associated with de- 
creased phosphorylation of the FGFR effectors extracellu- 
lar signal-regulated protein kinases 1 and 2 (ERKl/2) and 
with decreased staining for the proliferation marker Ki67 
(Figure 5E). Interestingly, both of the two most effective 
shRNA constructs elicited a compensatory increase in the 
expression of the a isoform of FGFR2 (Additional file 8: 
Figure S4). These results and previous data suggesting an 
antitumorigenic role for the FGFR-a isoform [15,16] led 
us to test the hypothesis that pulmonary metastatic out- 
growth is stimulated specifically by the expression of 
FGFR1-|3 and, conversely, is inhibited by that of FGFRl-a. 
To this end, we rescued FGFRl -depleted D2.A1 cells (that 
is, shFGFRl#2) with expression of a non-targeted form 
of FGFRl-a-IIIc (Figure 5F). When passaged at the 
same density as control cells, D2.A1 cells that expressed 
ectopic FGFRl-a displayed a decreased rate of cell 
growth (Additional file 10: Figure S6). 



FGFRl is known to be transactivated by the extracellu- 
lar matrix [38], and FGF signaling can be modulated by 
|33 integrin [39]. Along these lines, we previously dem- 
onstrated enhanced expression of |33 integrin in D2.A1 
cells when propagated within a three-dimensional extra- 
cellular matrix as compared to traditional two-dimensional 
tissue culture plastic [2]. Therefore, we sought to assess the 
impact of FGFRl-a expression on the three-dimensional 
outgrowth of the D2.A1 cells. Importantly, enhanced 
expression of FGFRl-a inhibited the three-dimensional 
outgrowth of D2.A1 cells to an extent similar to that 
achieved by depletion of total FGFRl (Figures 5F and 
5G). Consistent with these three-dimensional orga- 
notypic outgrowth analyses, D2.A1 cells that expressed 
ectopic FGFRl-a also demonstrated a dramatically re- 
duced capacity to produce pulmonary tumors in mice 
upon tail vein inoculation (Figures 6A and 6B), an 
event that significantly extended the lives of mice har- 
boring FGFRl-a tumors (Figure 6C). Collectively, 
these findings indicate that the antitumorigenic nature 
of FGFRl-a is dominant within the pulmonary micro- 
environment. Moreover, these data underscore the 
dichotomous nature of the a and |3 isoforms of FGFR, 
a relationship that can be understood diagnostically 
(Figure 4D) and used to guide the application of 
therapeutic inhibitors against FGFR for the treatment 
of BC. 




Figure 6 Expression of fibroblast growth factor receptor type ^, a isoform, has a dominant inhibitory affect during pulmonary tumor 
formation. (A) D2.A1 cells were constructed to express a scrambled short hairpin (shRNA; sc or scram) or the #2 shRNA targeting murine 
fibroblast growth factor receptor type 1 (FGFRl; sh or shFrl). Depletion of endogenous FGFRl in these cells was rescued with a nontargeted 
human full-length form of the third (III) extracellular immunoglobulin (lll-lg) domain of the a isoform of FGFRl (FGFRl -a-lllc; Frl -a) or green 
fluorescent protein (GFP) as a control. These cells (5 x 10^) were injected into the lateral tail vein of female BALB/c mice. Shown are representative 
bioluminescence images of mice from each group at the time of injection (TO) and 25 days later (day 25). (B) Bioluminescent quantification at the 
indicated time points for the cohorts described in (A). Data are the mean (+SE) thoracic area flux values normalized to the injected values, 
resulting in the indicated P values (n = 5 mice per group). (C) Survival analysis of the cohorts described in (A) resulting in at least the indicated 
P value (n = 5 mice per group). 
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Pharmacologic inhibition of FGFR kinase activity delays 
pulmonary outgrowth 

Given our findings thus far demonstrating the contrast- 
ing nature of FGFRl-a versus FGFRl-p in regulating 
pulmonary metastatic outgrowth, we next sought to as- 
sess the impact of therapeutically targeting the shared 
kinase domain of these FGFR isoforms. Importantly, 
FGF2-induced phosphorylation of ERKl/2 within tumor 
organoids growing in a three-dimensional matrix that 
recapitulates the pulmonary microenvironment could be 
effectively blocked using two different pharmacological 
inhibitors of FGFR kinase activity (Figure 7A). To model 
the therapeutic treatment of established pulmonary me- 
tastases, we grew D2.A1 cells for 4 days under three- 
dimensional organotypic culture conditions, at which 
point tumor organoids were treated with the clinically 
relevant FGFR kinase inhibitor, BGJ-398 [40]. Treatment 
with BGJ-398 significantly decreased the subsequent 
three-dimensional outgrowth of these established tumor 
organoids (Figure 7B). The importance of this finding is 



underscored by the fact that FAKl/2 kinases are known 
to be critically involved in the initiation of three-dimensional 
outgrowth [25,35]. However, treatment of these established 
and actively growing tumor organoids with the FAKl/2 in- 
hibitor PF271 (PF271) actually increased their outgrowth 
(Figure 7B). These findings, together with those of our previ- 
ous studies [2,25], highlight the critical changes that tal<e 
place as outgrowth initiation is driven by a mesenchymal 
state, which gives way to a MET program during macro- 
scopic pulmonary outgrowth. Importantly, the maintenance 
of MET-associated outgrowth requires FGFR signaling, but 
not FAKl/2 signaling. To verify the in vivo efficacy of 
BGJ-398 within the pulmonary microenvironment, D2.A1 
cells were injected into the lateral tail vein and allowed to es- 
tablish pulmonary lesions for 7 days, at which point mice 
bearing D2.A1 pulmonary tumors were treated daily with 
BGJ-398 via oral gavage. As shown in Figure 7C, administra- 
tion of BGJ-398 impaired pulmonary tumor outgrowth to an 
extent similar to that of PF271 (Figure 7C). Taken together, 
these findings highlight the complexities of outgrowth 




Days posi tail vein injection 



Figure 7 Therapeutic targeting of fibroblast growth factor receptor kinase activity delays the outgrowth of established pulmonary 
tumors. (A) Metastatic D2.A1 cells were grown under three-dimensional organotypic conditions for 4 days and not stimulated (NS) or pretreated 
with fibroblast growth factor receptor (FGFR) kinase inhibitor PDl 73074 (PD173) or BGJ398 (BGJ) for 6 hours prior to being not stimulated (ns) 
or stimulated with FGF2 for 30 minutes, at which point the cells were analyzed by immunoblotting for phosphorylated and total extracellular 
signal-regulated kinases 1 and 2 (pErkl/2 and tErkl/2, respectively). (B) D2.A1 tumor organoids were established for 4 days under three-dimensional 
culture conditions, at which point the tumor organoids were treated with BGJ398 (1 pM) or PF271 (1 |jM). Data are the bioluminescence values 
normalized to the plated value and are the mean (±SE) of two independent experiments completed in triplicate, resulting in the indicated 
lvalues. Inset: Photomicrographs (lOOx magnification) of the D2.A1 three-dimensional tumor organoids 4 days after initiation of treatment 
with the indicated inhibitors. (C) D2.A1 cells were injected into the lateral tail vein of female BALB/c mice and allowed to establish pulmonary 
tumors for 7 days, at which point the mice were randomized and split into the three cohorts that received vehicle, BGJ398 or PF271 
(50 mg/kg) daily via oral gavage (arrow indicates treatment initiation). Pulmonary tumor outgrowth was monitored by bioluminescence at the 
indicated time points (n = 5 mice per group). Bioluminescent images and ex vivo photomicrographs show the lungs from representative mice 
in each treatment group. Data are the mean (±SE) pulmonary bioluminescence values at the indicated time points relative to the injected 
value (% of DO), *P<0.05. 
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initiation and maintenance and suggest that therapeutic tar- 
geting of FGFR can potently inhibit these processes. 

Discussion 

EMT:MET cycles represent essential physiological pro- 
cesses that occur during critical points in the development, 
maintenance and repair of wounded epithelial tissues [41]. 
Through these processes, normal epithelial cells have the 
capacity to take on certain mesenchymal characteristics 
and then accurately return to their initial epithelial state. 
However, as has been observed in numerous other studies 
and quantified in our present study as shown in Figure 1, 
the EMT:MET process is pathologically engaged by carcin- 
oma cells during the metastatic cascade. Therefore, in the 
present study, we sought to address the hypothesis that, 
following initiation of metastatic outgrowth, tumor cells in- 
accurately complete the MET program and enter into a 
secondary epithelial state that is similar to, but critically 
unique to, the epithelial characteristics of the primary 
tumor cells from which they are derived. Using microarray 
analyses, we defined numerous factors that are differen- 
tially expressed between nonmetastatic mammary tumor 
cells and those that have undergone a metastasis-inducing 
treatment with TGF-p. Overall, these findings will contrib- 
ute to our understanding of the unique growth properties 
of metastatic lesions as compared to their corresponding 
primary tumors. 

Using the spontaneously metastatic 4T1 model in 
combination with a stable E-cad-luciferase reporter sys- 
tem, we quantified the dynamic regulation of E-cad ex- 
pression during the various steps of the metastatic 
cascade (Figure 1). Because these cells were derived from 
a spontaneous tumor, however, they almost certainly 
underwent one or more EMTiMET cycles during their 
original and natural development. Therefore, we utilized 
our EGFR transformation model to specifically demon- 
strate that TGF-|3 treatment is sufficient for the acquisition 
of metastatic properties (Figure 2). Indeed, immunohisto- 
chemical analyses and subculture of the primary tumors 
demonstrated that several changes in gene expression initi- 
ated in vitro were maintained in the primary tumor. How- 
ever, ex vivo subculture of the resulting metastases clearly 
indicates that these cells have renewed expression of E-cad 
and enter into a secondary epithelial state. Just as with the 
4T1 model, this new epithelial state easily gives way to a 
secondary spontaneous EMT and metastatic cycle, as sec- 
ondary engraftment of these cells onto the mammary fat 
pad leads to robust pulmonary metastasis without TGF-p 
treatment prior to their inoculation into mice. Subculture 
of these metastases yielded the NME-LM2 cell line and 
established an isogenic progression series of cell lines that 
possess increasing metastatic potential, ranging from nor- 
mal (NMuMG) to low-grade and nonmetastatic (NME) to 
high-grade and metastatic (NME-LM). 



Using our NME progression series, various human BC 
cell lines and patient-derived tissue samples, we verified 
the stable upregulation of FGFRl-p-IIIc that we initially 
identified in our in vitro TGF-p treatment and recovery 
microarray analyses. As such, we further investigated the 
role of FGFR in facilitating late-stage metastatic tumor 
outgrowth. Our data clearly indicate that cellular trans- 
formation is required in conjunction with TGF-p treatment 
to facilitate sustained upregulation of FGFRl (Figure 3). In 
this case, the means of cellular transformation is overex- 
pression of EGFR, and we previously demonstrated that 
EGFR is downregulated following the in vitro TGF-p treat- 
ment and recovery protocol used in our present study [2]. 
Therefore, it is interesting to note that maintained upregu- 
lation of FGFR may serve to explain the disparity between 
the power of EGFR as a predictive marker of poor progno- 
sis in BC [6] and the failure of EGFR- targeted therapies in 
the treatment of metastatic BC [7-9]. Studies aimed at fur- 
ther elucidating the role of FGFR in the inherent resistance 
of BCs to EGFR-targeted therapies are currently underway 
in our laboratory. Along these lines, we previously demon- 
strated that metastatic D2.A1 cells have diminished expres- 
sion of EGFR compared to their nonmetastatic and 
isogenic D2.0R counterparts [2]. Interestingly, D2.A1 cells 
predominantly express FGFRl-p-IIIc, which is consistent 
with their metastatic phenotype. Importantly, depletion of 
total FGFRl and ectopic expression of FGFRl-a-IIIc simi- 
larly inhibited pulmonary tumor outgrowth in the present 
study (Figure 6). Thus, our findings expand upon the 
results of previous studies that have linked FGFRl-p 
expression to the development of BC [15] and work in 
glioblastomas whose apoptosis was readily induced by 
administration of morpholino oligonucleotides to rees- 
tablish inclusion of the a exon [16]. 

Using flanking primer sets, we were able to identify 
human mammary tumor samples that not only upregu- 
lated FGFRl expression but also aberrantly excluded the 
a exon as compared to their matched normal samples 
(patient 1 and patient 7) (Figure 4). Developing this 
assay as a diagnostic screening test to detect individual 
FGFR isoform expression could prove to be highly bene- 
ficial in prospectively identifying those patients who 
would most likely benefit from FGFR inhibitor therapy. 
Indeed, administration of BGJ-398 completely inhibited 
the activity of FGFR under physiologic conditions and 
potently delayed pulmonary tumor outgrowth of the D2. 
Al cells. These data are consistent with the fact that the 
D2.A1 cells primarily express FGFRl- p. However, given 
our genetic studies identifying the antitumorigenic na- 
ture of FGFRl-a, administration of FGFR inhibitor ther- 
apies to patients whose tumors express this isoform 
(patients 2 and 13) (Figure 4) could potentially prove to 
be detrimental. In fact, this pro- and antitumorigenic di- 
chotomy between FGFR isoforms likely contributes to 
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Additional file 5: Figure S2. Transforming growth factor (3 treatment 
and metastatic subculture select for cells with an epithelial phenotype 
but enhanced metastatic potential. (A) Schematic representation and 
phase contrast photomicrographs of the normal mammary epithelial 
(NME) cell progression series. (B) Following transforming growth factor (3 
(TGF-[3)-induced metastasis, several markers of epithelial-mesenchymal 
transition (EMT) (epithelial cadherin (E-cad), vimentin, and neuronal 
cadherin (N-cadherin)) in the NME lung metastatic (NME-LM) cells 
returned to levels comparable to those of the parental NME cells. 
(C) RT-PCR analysis of E-cad showing similar levels across the NME 
progression series. Data in (B) and (C) are representative of two independent 
experiments yielding similar results. 

Additional file 6: Figure S3. Autocrine expression of epidermal growth 
factor receptor ligands are downregulated following transforming growth 
factor (3 treatment and recovery. Normal mammary epithelial (NME) cells 
were stimulated with transforming growth factor (31 (TGF-[3l; 5 ng/ml) 
and allowed to recover as detailed in the Methods section of the text. 
RNA was gathered, and global gene expression was assessed by 
microarray analysis (Supplementary Table SI). Real-time PGR was carried 
out to confirm downregulation of the epidermal growth factor receptor 
(EGFR) ligands betacellulin, amphiregulin, and epiregulin. Numbers indicate 
fold downregulation for each gene. 

Additional file 7: Table S4. Fibroblast growth factor receptor 
expression is increased in breast cancer tumors as compared to matched 
normal tissues. Twenty breast cancer (BG) tumor biopsies were analyzed 
by RT-PGR for expression of fibroblast growth factor receptor (FGFR) types 1 
to 4, whose expression levels were normalized against RNA gathered from 
adjacent normal tissues. Data shown are the mean (±SE) fold increases in 
FGFR levels as compared to the matched normal tissue. 

Additional file 8: Figure S4. Genetic depletion of fibroblast growth 
factor receptor type 1 abrogates pulmonary tumor formation in mice. 
Bioluminescent images of representative mice injected with D2.A1 cells 
expressing a control (scrambled short hairpin RNA; scram) or two unique 
fibroblast growth factor receptor type 1 (FGFRl)-targeting short hairpin 
RNAs (shFrl#2 and shFrl#3). Images were taken at the time of injection 
(time 0; TO) and 3 weeks following injection. 

Additional file 9: Figure S5. Genetic depletion of fibroblast growth 
factor receptor type 1 results in a compensatory increase in fibroblast 
growth factor receptor type 2, a isoform, expression. (A) Real-time PGR 
analysis of fibroblast growth factor receptor types 2 to 4 (FGFR2 to 
FGFR4) in D2.A1 cells expressing a nontargeting control scrambled short 
hairpin RNA (scram) or two unique short hairpin RNA (shRNA) sequences 
targeting FGFRl (shFGFRl#2 and shFGFRl#3). Data are the mean (±SD) of 
four independent experiments. (B) RT-PGR analyses of cells expressing 
FGFRl shRNAs as described in (A) using isoform-specific primer sets for 
FGFRl and FGFR2 as detailed in Figure 4A in the main text. 

Additional file 10: Figure S6. Expression of third extracellular 
immunoglobulin domain of fibroblast growth factor receptor type 1, 
isoform a, decreases the in vitro growth of D2.A1 cells. D2.A1 cells 
expressing a control scrambled short hairpin RNA (shRNA; scram) or a 
fibroblast growth factor receptor type 1 (FGFRl )-specific shRNA 
(shFGFRl) were transduced with retroviral particles harboring transcripts 
of either green fluorescent protein (GFP) or the third (III) extracellular 
immunoglobulin (lll-lg) domain of fibroblast growth factor receptor 
type 1, isoform a (FGFRl -a-lllc). These populations were selected with 
G418 and cultured at identical cell densities for an additional three 
passages (PI to P3). 



the limited efficacy of BGJ-398 in inhibiting the out- 
growth of the D2.A1 cells, as these cells do endogen- 
ously express detectable levels of the FGFRl -a isoform. 
In contrast to this cell culture model, the breast tumor 
samples of two patients analyzed in our present study 
(patients 1 and 13) yielded very clear upregulation of 
one isoform or the other, again supporting the use of 
this FGFR diagnostic approach as a predictive biomarker 
for initiation of FGFR inhibitor therapy. 

Conclusion 

Our studies demonstrate the dynamics of EMT and 
MET as BC progresses from carcinoma in situ to full- 
blown metastatic disease. FGFR and several other factors 
identified herein represent a signature of oncogenic 
TGF-|3 signaling that does not return to baseline during 
recovery from ligand exposure. This failure to accurately 
execute the MET process sets the stage for FGFR to act 
as a potent driver of pulmonary metastatic outgrowth, 
even if it may not have been an initiator of primary 
tumor formation. Overall, our findings have important 
implications related to the means by which science and 
medicine undertake targeting of FGFR for the treatment 
of metastatic BC. 

Additional files 



Additional file 1: Table 51. Table listing the sequences of the oligos 
used for the indicated applications. Also listed are the sequences of the 
fibroblast growth factor receptor type 1 (FGFRl )-targeted short hairpin 
RNA (shRNAs). 

Additional file 2: Table 52. Table listing genes that were regulated 
more than threefold between untreated normal mammary epithelial 
(NME) cells (before transforming growth factor treatment; pre-TGF) and 
those that had transitioned through a cycle of epithelial-mesenchymal 
transition (EMT) and mesenchymal-epithelial transition (MET) (EMT:MET) 
(post-TGF-Rec) as described in the Methods section of the text. As 
shown, microarray analyses were conducted in triplicate, resulting in 
the indicated averages and calculated fold changes in expression. The 
complete microarray data set can be found in the Gene Expression 
Omnibus database [GEO:GSE54491]. 

Additional file 3: Table 53. Table listing the antibodies used for the 
indicated applications, the dilution at which they were used and the 
supplier information. 

Additional file 4: Figure 51. Transforming growth factor (3 (TGF-(3)- 
induced normal mammary epithelial (NME) cell metastases are highly 
aggressive upon secondary fat pad inoculation. (A) Following primary 
tumor removal, untreated (Pre-TGF) and TGF-^-treated (Post-TGF) NME 
primary tumors were disassociated and subcultured in the presence of 
hygromycin (500 |jg/ml). The resultant cultures were grown on two-dimensional 
standard tissue culture plastic (2D-plastic) or under three-dimensional 
organotypic culture conditions (3D culture). Shown are representative 
phase contrast images (original magnification, 200x) depicting the 
typical growth morphologies of these ex vivo tumor cells. (B) Ex vivo 
tumors derived from TGF-(3-treated and untreated tumors were disassociated 
and analyzed by flow cytometry for expression of epithelial cadherin (E-cad). 
(C) Parental NME cells and their TGF-(3-induced lung metastatic derivatives 
(NME-LM; 1x10^ cells/mouse) were engrafted onto the mammary fat pad of 
nu/nu mice, and primary tumor growth and recurrence were monitored by 
digital caliper measurements (n = 5 mice per group). 
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